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I. INTRODUCTION 



One would like to have global solutions for rotating objects in general relativity consisting 
of a matter-filled interior region and a vacuum, asymptotically flat, exterior, with the aim of 
modelling planets, stars, star clusters, galactic nuclei or galaxies. The interior and exterior 
would be matched across a boundary, the surface of the object, E. Finding such global 
models is very difficult even for axially symmetric configurations in equilibrium. So far 
there are no explicit global models known other than those for spherical stars, which must 
be non-rotating, and Neugebauer and Meinel's disc of dust [1], which has no interior (the 
matter source has zero thickness and is described by jumps in the metric derivatives). 

There have been several recent studies of the structure of the underlying equations for 
the problem. The two regions can be treated independently subject to the matching at 
the boundary. For a given interior, this fixes "Cauchy data" giving both the metric and 
its derivative at the boundary, i.e. gives both Dirichlet and Neumann boundary conditions 
for the elliptic equations governing the vacuum exterior. A compatibility problem therefore 
arises. Uniqueness of the exterior solution has been shown in [2, 3] (independent of any non- 
circularity in the interior), and necessary conditions on the Cauchy data for the existence 
of the exterior also exist [4]. Correspondingly, the interior must satisfy these conditions 
in order to describe an isolated rotating compact object in equilibrium. Moreover, the 
conditions were found to be sufficient for static exteriors [5], and it has been conjectured the 
same holds for the stationary case, but this issue is still under investigation [4]. 

However, applying these conditions in specific situations has proved very difficult. For 
example, we had hoped to be able to use them to give a definitive treatment of the matching 
problem for the Wahlquist solution [6-8] but were unable to evaluate the required integrals 
in closed form. The understanding of rotating objects has advanced over the years through 
two major approaches, numerical relativity and perturbation theory. The latter has proved 
to be very useful and has been widely employed, leading for example to a well-developed 
theory of slowly rotating stars. Probably the key paper in the theory of slowly rotating stars 
is that of Hartle [9] (see also [10, 11]). He discussed the case of a rigidly rotating non-singular 
perfect fluid interior with reflection symmetry, with perturbations dependent only on the 
slow rotation imposed on a static background. The uniformity of rotation was justified by the 
argument that configurations minimizing the total mass-energy must rotate uniformly, and 
hence so must all stable configurations [9, 11]. He implicitly assumed spherical symmetry of 
the background (an assumption later shown to be true for many cases [12]), thus excluding 
convective motions in the interior, and the admissibility and C 2 character of the coordinates 
used, which is greater differentiability than required by the usual matchings [13]. 

With the aim of studying models for stars (not necessarily in equilibrium) using pertur- 
bation theory, there has been another direction of research. This has focused on the study 
of general perturbations around configurations constructed by the matching of static spher- 
ically symmetric spacetimes. The first works in this direction were the classical papers by 
Gerlach and Sengupta [14, 15], whose procedure was better justified by Martm-Garci'a and 
Gundlach [16] who considered gauge-independent quantities. The complete formulation of 
the relevant perturbed matching conditions to first order, and in general, was eventually 
given in independent papers by Battye and Carter [17] and Mukohyama [18]. Since these 
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works are not in principle aimed at isolated stars in equilibrium, there are no restrictions 
on the exterior due to asymptotic flatness or stationarity, and only first order perturbations 
are considered. The latter is important, since, as was pointed out in [9], we need to go to 
second order for isolated stars in equilibrium in order to obtain the effects of slow rotation 
on the shape of the star. The conditions to second order have been obtained only recently 
by one of us [19]. 

Yet another method used for the construction of models of slowly rotating stars has 
recently [20, 21] produced a particular model to second order in the approximation for 
a rigidly rotating and constant density perfect-fluid interior. The method is based on a 
two-parameter perturbation scheme, combining the post-Minkowskian and slow-rotation 
approximations, both taken up to the first non-linear level. The drawback of the model is 
inherent to the method, since there is no rotating exact Newtonian limit. 

In this paper we re-examine the bases of the perturbed matching theory. An obvious 
question is, therefore, what new things can be said at this fundamental level? In our opinion, 
there are several issues that need to be clarified. Firstly, have all the hypotheses involved 
in the analysis been spelt out in full detail? Secondly, if not, what are they, and are they 
indispensible? This is especially relevant in the matching of spacetimes, which lies at the 
very heart of any approach dealing with finite objects with boundary, such as models of 
stars. Thirdly, most, if not all, of the analyses in the literature deal with a specific matter 
model in the interior, usually a rigidly rotating perfect fluid, and the results do depend 
strongly on this assumption. It is natural to ask how one can develop a theory which is as 
independent as possible of the interior matter model (or even fully independent). Achieving 
this is clearly of interest, since the resulting theory could be applied to many different 
situations, from differentially rotating perfect fluids to convective fluids, viscous fluids, or 
even totally different matter models like elastic bodies. The aim of this paper is to answer 
all these questions in detail. 

One disadvantage of our approach, of course, is that without assuming a particular matter 
model we cannot study important questions like uniqueness or existence of interior solutions 
for given boundary data. However, even with specific interiors, the issue of existence of 
global models for self-gravitating isolated rotating objects is a very difficult one and ex- 
tremely few results are known. The only general theorem in relativity of this type to date 
refers to stationary and axially symmetric rigidly rotating fluids with spatially compact 
sections. Existence of such rotating stars was proven by Uwe Heilig [22] for relativistic con- 
figurations close to Newtonian models and sufficiently small (but finite) angular velocity. In 
perturbation theory the situation is not much better, although the results by Hartle [9] indi- 
cate a possible path towards a rigorous existence theorem of rigidly rotating perturbations 
of a fluid around a static spherically symmetric configuration. 

We believe that this paper is the first to give a comprehensive and consistent theory 
of rotating stationary and axisymmetric objects to first and second order in perturbation 
theory from first principles. We start from the Darmois form of the matching conditions, 
i.e. we do not a priori assume admissibility of particular coordinates, and specialize this 
form to the stationary axisymmetric case. We do not assume that the background about 
which we perturb is spherically symmetric (unlike previous work), though we shall give the 
specialization to this case, nor, as mentioned above, do we restrict the type of matter in 
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the interior region. The main restrictions are that: we consider only perturbations about 
static solutions; we assume that the exterior is an asymptotically flat vacuum and does not 
contain an ergoregion; and we assume that the axis is everywhere regular. 

It would in principle be possible to undertake a parallel study of the more general case 
where the initial configuration is stationary rather than static, but we have not explored this 
at any length. Since necessary conditions for the existence of a matching of a given interior 
to an asymptotically flat exterior can be found in the fully non-linear regime for this case, 
they can also be found for the perturbations. However, proving sufficiency, even if possible, 
appears to be considerably more difficult. 

In section II we give the general conditions for matching, specialized to the stationary 
axisymmetric case and expressed using Weyl coordinates in the exterior. In section III 
we introduce the perturbation scheme, in particular introducing some gauge choices. The 
perturbations in the exterior up to second order are considered in section IV: we give the 
perturbed Ernst equations together with the perturbed Cauchy data conditions resulting 
from the matching. Necessary and sufficient conditions on these data for a solution to the 
exterior equations are derived in section V (recall that sufficiency has not yet been shown 
for the corresponding full equations). These conditions are first expressed in terms of the 
equality of volume integrals over the exterior and surface integrals involving the Cauchy 
data and a principal function on the surface itself, the latter being defined as a solution to a 
partial differential equation (PDE). This form is then reduced to integrals on the boundary 
surface by first expressing the integrands as divergences of one-forms, which themselves are 
ultimately defined as solutions of ordinary differential equations (ODE). The final results are 
thus found in the form of integrals over the parameter which runs along the meridians of the 
boundary surface. Lastly we specialize the static background to the spherically symmetric 
case in section VI: here some of the functions can be found explicitly, since the exterior 
background must be the Schwarzschild solution. We note that we do not in this paper use 
the assumptions of Hartle's work [9, 10]: we intend in a later paper to spell out fully the 
relationship of his approach and ours. 



II. MATCHING WITH WEYL COORDINATES 

In our problem the spacetime is composed of two regions with boundary, namely the 
interior (V 7 , g 1 ) and the exterior (V E , g E ) (where I and E stand for interior and exterior 
from now on). Each region, as well as its boundary, is assumed to be stationary and 
axially symmetric, and to have a regular axis of symmetry (see e.g [23]). We make no 
specific assumption on the energy-momentum tensor of (V 1 , g 1 ): not even the so-called 
circularity condition [24] in the interior will be assumed (which for fluids without energy 
flux would be equivalent to the absence of convective motions [25, 26]). In other words we 
shall not assume that the orbits of the stationary axisymmetric isometry group in the interior 
region are orthogonally transitive (i.e. the metric might not admit coordinates adapted to 
the isometries in which the line-element becomes block diagonal). We do not, however, 
treat interiors containing non-gravitational fields which propagate to the exterior (such as 
electromagnetic fields, for instance). This means that the exterior region will be taken to be 
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vacuum. Moreover, all matter configurations will be taken to be isolated, and the spacetimes 
are assumed to be asymptotically flat. 

Let us denote by X/ and E E the respective boundaries of (V 1 , g 1 ) and (V E , g E ). For 
us to be able to identify them, these hypersurfaces must be diffeomorphic to one another. 
As usual, this is best described by taking an abstract three-dimensional manifold a and 
two embeddings x 1 and x E ', where x I/E ■ o -»• such that = x I/E (°)- The 

point-to-point identification of the boundaries is x 1 ° (x E )~ 1 defined on S E . 

We assume also that in general the spacetimes admit no further local symmetries beyond 
stationarity and axial symmetry (though when perturbing we shall allow the background 
spacetime to be spherically symmetric). Then the axial Killing vector is uniquely fixed by 
demanding that its orbits are closed and the Killing coordinate has periodicity 2n. Let 
us denote by ff 1 and ff E these unique axial Killing vectors in the interior and the exterior 
spacetimes, respectively. As (V E , g E ) is asymptotically flat, it admits a unique Killing 
vector, denoted by C, E , which is unit timelike at infinity. In the interior region there is 
no equivalent way of fixing the stationary Killing vector uniquely before the matching is 
performed. Two different situations might be considered, one in which the interior (V 1 , g 1 ) 
is unknown and needs to be determined, and the other where it is explicitly known. In 
either case, we can choose £ to be the unique Killing vector which matches continuously 
with £ E (recall that the boundaries are also assumed to be stationary and axially symmetric 
so that we have a "symmetry-preserving matching" [27]), i.e. we can propagate the exterior 
Killing vector into the interior, but in the second case this need not agree with the timelike 
Killing vector used in writing g 1 , so we have to introduce two extra parameters in the interior 
metric in order to allow the most general matching. This is well understood and will not 
be discussed further here; see [2]. The isometry group being two-dimensional and the orbits 
of fj 1 ^ closed, the Killing vectors £ J / E and fj 1 ^ commute ([28], see also [29]). Thus, away 
from the symmetry axis, there exist local coordinates {T, <3>, xf} (A, B ■ • ■ = 2, 3) in (V 1 , g 1 ) 
such that I; 1 = dx and fj 1 = <9$. 

For the exterior metric, we do not allow ergoregions (as we aim to model stars rather 
than black holes) and the vacuum field equations then imply the local existence of Weyl co- 
ordinates {t, 0, p, z} which satisfy the conditions that (i) they are adapted to the isometries, 
i.e. £ B = d t and fj E = 8$, (ii) p = defines the axis of symmetry and (iii) the metric g E 
takes the local form 

ds| = -e 2U (dt + Ad<j)) 2 (1) 
+e~ 2U [e 2k (dp 2 + dz 2 ) + p 2 d0 2 ] , 

where U, A and k are functions of p and z only. The scalar p is intrinsically defined as 
the non-negative solution of p 2 = — (£ £ , (, E ) g E(fj E , fj E ) g E + (£ £ , fJ E ) 2 gE where (,) g denotes 
the scalar product with respect to a metric g. z is also intrinsically defined as the scalar 
(unique up to a sign and a constant shift) such that dz is orthogonal to dp (with respect 
to the metric g E ) and dz and dp have the same norm. The coordinate freedom in (1) 
consists only of constant shifts of t, and z. For later use, let us recall that stationary and 
axially symmetric vacuum spacetimes admit, locally, a scalar function Q (the twist potential) 
defined, up to an additive constant, by 

P n p = -e w A, z , P n z = e AU A p . (2) 
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If, moreover, (V E , g E ) is simply connected Q can be defined globally. In asymptotically 
flat spacetimes the additive constant is fixed so that Q — > at spatial infinity. The scalars 
U and Q are intrinsically defined by e 2U = —{!; E ,!; E ) and df2 = *(£ E A d£ E ), where 
£ E — g E (C i E ,-) and * denotes the Hodge dual [39]. The Einstein vacuum field equations 
reduce to a complex second order PDE for (U, Q) (a pair of real PDEs in the form used 
below), the Ernst equation(s), together with quadratures for the metric functions A and k. 

With this choice of coordinates in (V E , g E ), and since the stationary and axial Killing 
vectors are tangent to the boundary X 7 , it follows that there exist local coordinates {r, ip, //} 
on the abstract manifold a and functions T ([/,), and xj(fi) (see [3], and [2] for the 

orthogonally transitive case), such that the embedding x 1 reads 

X 1 :{T = r + T(/i), $ = <p + $(//), xj = 

and, moreover, dx s (<9 T ) = £ 7 |e/, dx £ (<9^) = 77 7 1 S -r and dx E (d fJi ) is orthogonal to £ 7 and 
?/ 7 on X 7 . It can easily be checked that {r, </?,//} is defined uniquely except for constant 
shifts of r and y? and redefinitions fi(fi'). In terms of these coordinates on a, the exterior 
embedding x E is forced by the matching conditions (more precisely, by the continuity of the 
first fundamental form) to take the following form 

X E ■ {t = T, <f) = <p, p = p(/i), z = z(n)} 

for some functions and z(fi). Constant shifts in t and are in principle allowed, but 
they can be set to zero without loss of generality by exploiting the freedom in the exterior 
coordinates. With the embeddings at hand, we can discuss the necessary and sufficient 
conditions for the interior and exterior spacetimes to be joinable (i.e. to produce a spacetime 
(V, g) with continuous metric and no surface layers in the Riemann tensor [13]). These are 
the well-known Darmois conditions [30] which require that the first and second fundamental 
forms inherited on a from both sides agree. In our case, this set of matching conditions 
can be conveniently rewritten as follows (see [31, 32], [2] for the case where the interior is 
orthogonally transitive, and [3] for the generalization to an arbitrary interior): 

(a) Conditions on the interior hypersurface, given by the Israel conditions [33] 

n Ia n^S^ y = 0, n Ia e/ P S^ | s , = 0, (3) 
2 = 1,2,3, 
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where S 1 ^ is the Einstein tensor of g 1 , n 1 is a normal vector to X 7 in (V 7 , g 1 ) and e. 
are any three independent vectors tangent to X 7 . In principle, any choice of normal 
vector n 1 is suitable for writing (3). However, for the remaining matching conditions 
it is convenient to fix this vector so that it points to the interior of (V 7 , g 1 ) and has 
norm 

(™'."'),i,= (^e') s 4,, (4) 

where e 7 = dx / (9 M ). Note that redefining fj, on a changes e 7 and hence n 1 . Neverthe- 
less, all expressions below are easily checked to be invariant under rescalings //(//)• 
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Conditions (3) determine which (if any) hypersurfaces in a given interior metric are 
candidates to match to an empty exterior. If the matter content inside is, for instance, 
a perfect fluid with no convective motion, these conditions reduce to p — 0, where p 
is the pressure of the fluid. Generically (3) consists of four independent conditions, 
which means that only for special interior metrics will hypersurfaces matching with 
vacuum exist. Whenever the circularity condition is satisfied in the interior, two of 
(3) become trivial. 

(b) Definition of the exterior matching hypersurface. The functions p(p) and z(p) deter- 
mining the form of the exterior surface are uniquely fixed by 

p(p) = a kh z(p) = -r^(a)| s /, (5) 

where a 2 = — (£ 7 , £ 7 )^ (f) 1 , i7 7 ) s ' + (£ 7 , V 1 )^^ « > 0, and the dot denotes the derivative 
with respect to p. The additive constant in z(p) is inessential given the shift freedom 
z — > z + const. 

(c) Boundary conditions for the exterior problem. The rest of the matching conditions 
provide the following data on the exterior metric functions U and A on S s 

Uy = Vy, n E (U)y=n I (V)y, (6) 
A\vb = B\v, n E (A)y=n I (B)\ J]I} (7) 

where n E = -zdp+pd^E, e 2V = -(£ 7 ,£ 7 ) s /, and B = (£ 7 , ?7 7 ) g //(f 7 , f 7 ) 9 /. In order 
to ensure that n E points to the exterior of (V E , g E ) we choose p to be an increasing 
function of p at the south pole (i.e. at the intersection of S 75 with the symmetry 
axis having minimum value of z, which we assume exists). This choice constrains the 
allowed rescalings p(p') to be strictly increasing. 

The conditions for A translate into boundary conditions for the twist potential as 
follows 

0| E E = Vt^pP + Q, z z\y,e 



p4t/ 

— n E {A) 
P 

-n 7 (5) 



a 

3 4C7 
P 



a 



(8) 



EJ dp 



The right hand sides of (6), (7) and (8) are known once the interior is known. Thus 
the matching conditions fix the normal derivative of Q on S 75 uniquely and Q on H E 
up to an additive constant. This constant is in principle relevant as Q has been defined 
so that it vanishes at infinity. However, it can be proven [2] that there is at most one 
value of the additive constant for which the exterior vacuum field equations with the 
(overdetermined) boundary data (6) and (8) are compatible. 
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III. THE PERTURBED MATCHING CONDITIONS 



As usual in perturbation theory, we consider a one-parameter family (V t ,g t ) of 4- 
dimensional spacetimes, differentiable in e, and think of perturbations in terms of derivatives 
of the metric with respect to e, evaluated at e = 0, which requires working on a single man- 
ifold (more precisely, one manifold for the interior region and one manifold for the exterior 
region). Quantities in the previous section, except for the coordinates, as we shall discuss 
below, will now bear a subscript e. 

Let us discuss the interior region (the exterior case will be similar). In order to deal 
with a single manifold, we need to identify in some way points of different spacetimes V/ in 
the e-family. If the manifolds were without boundary and diffeomorphic to each other we 
could take any diffeomorphism, smooth in e, between, say, Vq and V/ in order to identify 
them. It is clear that such an identification is not unique, and that there is no canonical 
choice, because we can perform a diffeomorphism S e of V{ onto itself, before applying the 
diffeomorphism above. This freedom in performing the identification is the heart of the 
gauge freedom inherent to perturbation theory (see e.g. [34]). 

Once an identification has been chosen, we have a single manifold Vq and a collection of 
metrics g\ defined on it. In perturbation theory (up to n-th order), only the background 
metric g\ and the first n derivatives of g\ evaluated at e = are of interest. Thus we have a 
background spacetime (Vq, g^) and n symmetric tensor fields, (a = 1, • • -n), defined on 
it (the perturbations). The freedom in the identification of spacetimes translates into the 
gauge freedom in perturbation theory and is defined by n vector fields on Vq. For instance, 
to first order, and denoting by K[ the first order perturbation tensor, the gauge freedom is 
K'l = K[ + £g ± go where s\ is the first order gauge vector and £ denotes the Lie derivative. 
For higher order perturbations, the gauge transformations are more complicated as they also 
involve all lower order terms. 

These issues are all well-understood and would have deserved no inclusion here except 
that in our case the manifolds we are considering are with boundary. How to define the 
identification in this case has recently been discussed in [19]. We repeat the main idea here 
for completeness. For each e, V/ is a manifold with boundary E^. Thus identifying them 
via diffeomorphisms requires, strictly speaking, that boundaries are mapped into bound- 
aries. However, if we view each manifold V/ as a closed subset of a larger manifold without 
boundary, the condition that the boundaries are mapped to each other strongly restricts the 
gauge freedom (at least near the boundaries) and this may not be suitable for the problem 
at hand. It is more convenient to let the boundaries "move" freely in the identification. 
Perturbation tensors can still be defined everywhere on the background spacetime Vq with 
boundary Eq as follows. For points away from the background boundary Eq the usual pro- 
cedure obviously works. So we only need to worry about how to define perturbations on Eq. 
At each point p G Eg and for small enough positive e, the corresponding identified point in 
V/ will lie on E^ or move towards the interior of V/ or move towards its exterior (within the 
larger manifold). In the first two cases we can define the perturbation tensors by defining 
the derivatives as one-sided limits with e — > 0, e > 0. In the third case, the point p moves 
towards the interior of V/ for negative values of e. So derivatives can again be defined if we 
take the one-sided limits with e < 0. Note that this construction is independent of the larger 
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manifold with boundary, which can be dispensed with altogether. Hence we can, as before, 
define n symmetric perturbation tensors K T a on the background manifold with boundary 
(Vq, c/q) up to and including the boundary. Obviously, exactly the same considerations hold 
for the exterior region. 

Following the discussion above, we can introduce coordinates {T, &,xf} for each space- 
time (V/, g{). In principle, the coordinates themselves should have an e label. However, we 
can make this unnecessary by propagating the coordinates from (Vo, go), using the identifi- 
cations, so that points with the same coordinate values {T, Q,xf} in the interior regions of 
different spacetimes are identified. This choice reduces the gauge freedom available since we 
are only left with the changes of identification describable by coordinate changes of the type 
{T —>T + Ti{xf, e), $ — >■ $ + $i(xf, e), xf — > e)}. This (partial) gauge fixing is in 

fact very useful in the sequel. In this gauge, let us denote by K[ and K{ the first and second 
perturbations of the metric tensor on the interior background (Vq, g^), with boundary Eq. 

In the exterior regions, we analogously relate the spacetimes (Vf, gf) so that points with 
the same Weyl coordinates {£, 0, p, z} in different spacetimes are identified. Since the Weyl 
coordinates are unique except for constant shifts in t, <p and z, the freedom in performing 
the identification is reduced to t — > t + (3o(e), <j> — > <t> + Pi(c), z — > z + (3 2 (e). Thus the gauge 
freedom in the exterior region is reduced even more strongly than in the interior, which is 
also useful. The exterior background is therefore (V(f, go) with metric 

go = -e 2U °dt 2 + e~ 2U » [e 2k » (dp 2 + dz 2 ) + p 2 d0 2 ] , (9) 

and the perturbation tensors (up to second order) take the following form 



Kf = -2e 2U «U' dt 2 - 2e 2Uo A' dtd<P + 2e" 2l V fc ° {-% + k' ) (dp 2 + dz 2 ) - 2 e - 2C/ °^p 2 d0 2 , 



K 



-2e 



2U Ijjn 



2U' Q 2 ) dt 2 - 2e 2Uo (Aq + AA'M dtd0 - 2 



; 2C/ M' : 



(10) 



+ e" 2 V 



fa - 2^ 2 )] d<p 2 + 2e~ 2U °e 2k0 [k% -U^ + 2(k' - %f \ (dp 2 + dz 2 ) , 



which follow directly from (1) by taking e derivatives. Since a gauge choice is involved in 
the definition, we shall call them perturbation tensors in Weyl form or in Weyl gauge. In 
(9) and (10), Uq, k , U' Q , A' , k' , Uq, A' ' and k' ' are functions of p and z and are defined by 



Uq — U e \ t= o, Uq 



d e U e \ e -_ 



=0) u o 



d 2 e U e \e- 



=o, etc. 



IV. THE SECOND ORDER PERTURBED EXTERIOR PROBLEM 

Let us start by recalling the vacuum field equations for (1). As already mentioned, we 
only need to concentrate on the equations for U t and Q e , since the remaining field equations 
for k t and A t reduce to quadratures. 

Let us consider Euclidean space (E 3 , 7) in cylindrical coordinates, so 7 = dp 2 +d,2 2 +p 2 d(/> 2 . 
For each value of e, let us consider an axially symmetric surface S e defined by Xe = {p = 
p e (p), z = z e (/i), = <p} with the same functions p e (p) and z e (p) as those defining Sf in Vf . 
Let us denote by D € C E 3 the exterior region of S e . The use of Weyl coordinates in (Vf, gf) 
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allows us to consider {U e , Q e } as fields on D t . The vacuum field equations are equivalent to 
the so-called Ernst equations on D t 

A 7 f/ e + i e - 4 ^(da,da) 7 = 0, 
A 7 a - 4 (da, dC/ e ) 7 = 0, 

where A g is the Laplacian of the metric g. This set of equations is supplemented with the 
asymptotic values U e — > 1, a — > at infinity plus the boundary data on S e coming from 
the matching conditions. As already mentioned, these data determine both {U e , Q e } on E € 
and its normal derivatives (except for an additive constant in Q € ). Thus we are dealing with 
Cauchy data for the elliptic system of equations (11). This is an overdetermined problem 
and we should not expect solutions to exist for arbitrary data. That expresses the fact that 
given an arbitrary stationary and axially symmetric interior metric (even if it is perfect fluid, 
say), there will in general be no stationary and axially symmetric vacuum exterior solution 
matching with it and also asymptotically flat. Thus existence for the exterior problem is an 
important issue. This being true for all e, the same will happen for the perturbed matching 
and field equations, as we discuss next. In the following subsections, we first derive the 
perturbed field equations (subsection IV A) and then the boundary conditions (subsection 
IV B). Section V is devoted to discussing under what conditions the Cauchy data of the 
perturbed field equations are compatible. 



A. The Ernst equations up to second order 

Let us obtain the systems of equations satisfied by the different orders in e of U t and Q t . 
Differentiating (11) with respect to e we get 



A 7 C/ e ' + e" 4 ^ (da, dfi'J 7 - 2e- AU *U' e (da, dtt 

A 7 a e - 4 (da e , dc/ e ) 7 - 4 (da, du' £ \ = o. 

Differentiating once more we obtain 



0, 



(12) 



A 7 C/ e " + e 



-AUe 



su' e (da, dn f e )^ + (da, tax + ( dn 'e> dn 'e)- 



-2u" (da, da) 7 + 8tr (da, da) 7 
A 7 a; - s (da e , du'x - 4 (da;, dux - 4 (da, du'j) 



0, 
0. 



(13) 



The systems for the zeroth, first and second order are now obtained by evaluating the 
systems (11), (12) and (13) at e = 0. As already mentioned, we are interested in studying 
perturbations of static objects. Thus the background exterior metric satisfies a — and 
hence Uq is a solution of the Laplace equation 



AM = 0. 



(14) 
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In fact we are primarily interested in spherical backgrounds so that the exterior background 
metric is the Schwarzschild metric and Uq is the corresponding Schwarzschild solution. Nev- 
ertheless, for the sake of generality we shall keep an arbitrary static and axially symmetric 
background until explicitly stated. The equations for the first-order perturbation of the 
static background {Uq, £1' }, which will be called first-order linearized Ernst equations, read 

A,n;-4(dsj;,df/„) 7 = o, 

while for the second order perturbation {Uq, SIq} we get the second order perturbed Ernst 
equations 

A 7 Qq - 4 (dnl dlU - 8 (d%, dU' Q \ = 0. 



B. Cauchy boundary data up to second order 



In order to obtain the boundary data for the exterior problem we should, in a general 
setting, consider two background spacetimes (Vq, gl) and (Vq, g E ) which match across the 
unperturbed boundaries Eq and Ef according to the standard matching conditions. More- 
over, the first and second order perturbed metric tensors K T a and (a = 1, 2) should satisfy 
suitable conditions on Eq and E^ coming from suitable first and second e-derivatives of the 
full matching conditions. These conditions may be called first and second order perturbed 
matching conditions. Their explicit form in full generality has been obtained by one of us 
[19]. These can be specialized to find the perturbed matching conditions of interest in this 
paper. However, for the sake of self-consistency we shall follow an alternative procedure 
which is well adapted to the stationary and axially symmetric problem we have at hand. 
Since, as we saw in Sect. Ill, the boundary conditions for all e reorganize themselves into an 
elegant form and the exterior problem reduces to the Ernst equation for (U e , Q e ) alone, we 
need only concentrate on the perturbed boundary conditions for these objects. Note that we 
have introduced coordinates {r, ip, ji}, with no e label, on the abstract matching manifold 
a e . The reason is the same as before, i.e. we identify points in different a e having the same 
coordinate values {r, <p, ji}. 

Given the interior family of metrics, the matching conditions give us (for every e) two 
functions p e (/u) and z e (/i). Thus we have a family of axially symmetric surfaces E e in E 3 , all 
of them diffeomorphic to a . To avoid repetition, we collect together here the assumptions 
and notation we shall use. 

Assumptions: Let (Vq, (?o) an d (Vq , g E ) be two static and axially symmetric space- 
times which can be joined across their static and axially symmetric boundaries Eq 7 E^. 
Let (V(f, g$) be vacuum and asymptotically flat and choose, locally, Weyl coordinates (9). 
Let K 1 ^ and K^ E be symmetric tensors on V$ E , invariant under the static and axial 
isometries and such that K^ 2 take the Weyl form (10). Take any metric g\ in Vq such 

that g\ — g- + eK[ + \^K\ + 0(e 3 ) and denote by and fj^ its stationary and axial 
Killing vectors. Let (V/, g{) admit a hypersurface E^ where it can be locally matched to a 
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vacuum exterior. Let the identification of with (T be made through the common use 
of coordinates {p, </?} in both manifolds. On a let functions Pq(p), Pq{p), z' (p), z'q(p) be 
defined by p € (p) = p + ep' + |e 2 p ' + 0(e 3 ) = a e | s / ; z e (p) = z + ez' + \e 2 z'^ + 0(e 3 ) = 
— nf(a e )| S i, {where a e and n/ are as in (4) and (5), evaluated at S e ) 7 and functions V t , nV e , 
W e , nW e ,V (p), VM, V^'(p), nV (p), nV^(p), nV>>(p), W (p), W>(p), W>>(p), nW (p), 
nW^p), nW>'(p) by 

K = K. + e\/ ' + i e V " + O(e 3 ) 
nV t = nV + enV^ + -e 2 nV£ + 0(e 3 ) 

W e = eW> + ^e 2 WZ + 0(e 3 ) 

nW £ = enWo + -e 2 nW^ + 0(e 3 ) 

Note that the staticity implies that the background values are Q = W (p) = and 
nQo = nWo{p) = 0. As shown in the previous section, the matching conditions imply that 
V e , nV e , W e and nW e are provided by the interior metrics, and give the boundary values 
^e|s e = K, rt e (U e )\s e = nV e , fi e | Se = W e and n e (f2 e )|s e = nW t for the exterior problem. 
Note that the right hand sides of the four equations just given are functions of p and e 
alone, while the left hand sides are functions on E 3 evaluated on a (moving) surface £ e . We 
can now take derivatives of these expressions with respect to e (at constant p) in order to 
obtain the perturbed boundary conditions for the perturbation functions U' 0) Q' , U{{ and Qq. 
It is clear that when taking derivatives on the right hand sides of (17) two types of terms 
appear, namely those coming from the explicit dependence on e in the functions and those 
coming from the fact that the surfaces S e also depend on e. The latter will involve d p and 
d z (and higher) derivatives of lower order terms U , Q , Uq and Q' , all of them evaluated 
at the unperturbed surface So- Our aim is to express everything in terms of the twelve 
functions V (p), V^), V^'(p), nV (p), nV>(p), nV£(p), W (p), W>(p), W>'(p), nW (p), 
uWq(p), and nW^p) of p, which correspond exactly to the boundary information coming 
from the interior through the background matching and the first and second order perturbed 
matching conditions. 

Let us take a fixed point p e <xo and consider the trajectory in E 3 defined by Xe{p) when e 
varies. The first and second derivatives along this trajectory, evaluated at e = 0, give us two 
vectors on p and hence two vector fields on So- They determine how the matching surface 
moves within E 3 to second order in approximation theory. Denoting these vector fields by 

— * — * 

Z\ and Z 2 , we obviously have 

Z x = p' d p + z' d z \ So , Z 2 = pld p + 4U|e , (18) 
where, as before, p' = d e p t \ e=0 , p ' = d t d t p e \ t=0 , etc. Two further relevant vectors on S are 



2(&,£)rf s / 

d ^\{ii,ii) g i) 
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the first and second perturbations of the normals n e . Explicitly 



dn e 



de 
d 2 n f 



de 2 



6 = 



6 = 



-%)dp + p' d g \v , 



(19) 



Since we want to rewrite everything in terms of intrinsic objects on the unperturbed surface 
Sq, we need to use a basis adapted to this surface. A convenient choice consists of the 
unperturbed tangent and normal vectors, namely e = dxo(d^) = pod p + io<9z|s and n = 
n t \ e=0 (for notational convenience we drop a subindex both in e and in n). Writing down 
the vectors {d p , d z } on E in terms of e and n we get 



0. 



Z En 



d, 



z o + Po 



z 2 o+pl 



[z e + p nj , 
(^p e - z nj . 



(20) 



We can now express Zi, Z 2 , ^| e=0 , , in terms of {e, ft}. When doing this, six scalar 

fields on S appear in a natural way. Denoting them by Pi, Qi, P 2 , Q 2 , X and X 1 we have 



Z 1 = Pie + Qin, 



dn f 



de 
d 2 n e 



de 2 



6=0 



6=0 



where 



Pi = 



PoPo + zqZq 
Po + z o 



Z 2 = P 2 e + Q 2 n, 



dQ 



dp 
dQ 2 



dp 



Qi = 



- + P 1 X -Q 1 X 1 j e + 
+ P 2 X -Q 2 X 1 )e + 



dP 1 
dp 

dp 



+ PiX x + Q 1 X j n, 
+ P 2 X X + Q 2 X ) n, 



Po z o — z oPq 

• 2 i -2 ' 

po + 4 

PqZq — ZQpQ 

• 2 i • 2 : 

Po + z o 



p-l = 



X l 



PoPo + Z Z Q 

po + 4 

PopO + ZqZq 

pI + z I 



Q 2 = 



Po z o — z oPo 



Pi + z l 



(21) 



Having defined these objects we can evaluate the first and second order perturbations of the 
matching conditions. The statement of the result is rather lengthy. 

Proposition IV. 1 Under the Assumptions stated above, the metrics g{ and gf = g§ + 

eKf + \e 2 K2 match perturbatively to second order on T,^ E if and only if the following 
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conditions are satisfied 



U 



olsf 



n(Uo)\ 



^0 



Q' 



^0 



V , n(U Q )\ 
nV> + ^(q 



yE 



nV , 



U' 



olsf 



V> - P^ - Q 1 nV , 



dV \ d(Pi(n^ )) 



dp 
d 



+ Qi 



= Vq" -2P 1 -±- 2Q 1 nV ' + — (P 2 - Q\) 



dp 



dp 



po dV z 

— nV 

Po dp po 

dV \ d(2P 1 Q 1 nV ) 



dp 



dp 



f ( - P 2 f /' 2 .V: + 2AQ 1 X - Ot'.V, - Q 2 ^) ^ 

Po/ «P 

+ ( ~Qi - PiXo + 2P 1 Q 1 X 1 + g 2 x + Q\j^ ) n\/ , 

2 ^gK) , 2Qi fpodV' 



d 2 



dp 2 
d 

dp 
d 

dp 



2P 1 Q 1 



dp 
dV 



d 2 



dp J dp 2 



dp 

(P 2 -Ql)nV c 



Po dp po 



Q 2 + {P 2 - Ql) X - 2P 1 qA - 2P 1 Q 1 X 1 - Q{ 

Po Po 



2 Z 



dVo 
dp 



-Pi + {Pi - Ql) X 1 + 2P l Q l ^ + 2P 1 Q 1 X - Ql— 

Po Po, 



V[Q 2 + {Pi-QI)Xo-2P 1 qA 

Po 



f I 2P 1 Q 1 X - Q 2 ^ ) ( ^nV 4 



Po 



Po' 



£o 
Po 

Po dV 
Po dp 
zp dVp 
Po dp 



nVo 



—nVo 
Po 



dW' 



nWp + 2 



+2Q X 



- 2 



dp 



dp 

dVo\ d\Y!, 



o 



dp 
Po 



+ AnVo HWq 



where Pi, P 2 , Qi, Q 2 , X ± , X 2 are defined in (21). 

The actual necessity and sufficiency arises from the corresponding properties of the gen- 
eral Darmois conditions. The main work in the proof is the direct but cumbersome calcula- 
tion needed to arrive at the formulae: we leave this to the Appendix. 



V. COMPATIBILITY CONDITIONS 

In the previous section we found the overdetermined boundary data for Uq, Q' , Uq, Q' ' 
in terms of the interior metric and perturbation tensors. Since the equations they satisfy 
are elliptic we need to determine under what conditions asymptotically flat solutions exist. 
Asymptotic flatness requires that, for any value of e, lim p 2 +2 2^ 00 U e — 1 and lim p 2 +z 2_ >00 Q t = 
0. This implies lim p 2 +z 2_ >00 Uq = lim p 2 +z 2^ 00 Q' Q = lim p 2 +2 2^ 00 Uq = lim p 2 +z 2_ >00 Q' Q ' = 0. 
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The problem of determining the Neumann data (i.e. the normal derivative of the func- 
tion on the boundary dQ of a domain Q) in terms of the Dirichlet data (the value of the 
function on dQ), so that the Cauchy problem for an elliptic equation on Q is solvable, is 
called obtaining the Dirichlet-Neumann map. Many results are known on this problem, 
including deriving general properties of the map, determining the coefficients of the elliptic 
equation from the Dirichlet-Neumann map and finding explicit representations of this map 
for especially simple equations and domains. Good references on this topic are [35] and [36] . 
Our problem can then be phrased as saying that we want to find explicit restrictions on the 
boundary data so that they solve the Dirichlet-Neumann problem for the linearized Ernst 
equations on axially symmetric domains. 

While the equation for Uq is just the Laplace equation in Euclidean space, the equation 
for Q' has, in addition, lower order terms. In order to treat all cases at the same time it 
turns out to be convenient to define a conformally flat metric 7 = e~ 8C7 °7 on Dq, in terms 
of which the second equations in (15)-(16) can be rewritten as 

A^T = 0, 

A^' ' = 8(dfi' ,d^)~. 

Thus all equations for Uq, Uq, fl' Q , Q' ' can be collectively written as 

A^u = j, (22) 

where u = u(p,z) stands for Uq, Uq, etc., and j = j(p,z) represents the inhomogeneous 
terms in the second order perturbation equations. The metric 7 corresponds to either 
7, for the U -equations, or 7, for the ^-equations. The domain (Do, 7) is clearly non- 
compact because 7 is an asymptotically flat metric. Thus the compatibility conditions for 
the boundary values of Uq, Uq, £1' , Qq can be studied as particular cases of the compatibility 
conditions of the Cauchy problem for the general inhomogeneous Poisson equation (22) 
defined on a non-compact asymptotically flat region (Dq, 7) with a boundary dD = {p = 
Po(fj), z = zq(p), 4> = <p}. Furthermore, we shall assume that j tends to zero at infinity at 
least like 1/r 4 where r = a/ p 2 + z 2 (this requirement is fulfilled in the cases we are concerned 
with due to asymptotic flatness). 

We start with some well-known facts from potential theory. A simple consequence of 
Gauss' theorem is Green's identity which, for any compact domain K C Dq with C 1 bound- 
ary dK and any C 2 function ip on K with C 1 extension to K U dK, reads 

(t/jA^u — uA^t/j) rjj — ipn^(u) — un^(ip) dSj, 

K JdK L 

where is a unit (with respect to 7) normal vector pointing out from K, rj^ is the volume 
form of (D , 7) and dS^ is the induced surface element of dK. We intend to apply this 
identity to a function ip that (i) solves the Laplace equation A^ip = on Dq, (ii) admits a 
C 1 extension to dD and (iii) decays at infinity in such a way that ip^Jp 2 + z 2 is a bounded 
function on Dq. A function ip satisfying these three properties is called a regular ^-harmonic 
function on D (if a function satisfies just (ii) and (iii) and is C 2 on Dq we shall call it regular). 
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For such a function we can take K = D in (23) because the integral over the boundary "at 
infinity" can be easily shown to vanish. Thus 



ipn^(u) — ufij(ip) dSj = / ipjrjj. (23) 
*dD L J J D 

This expression relates the overdetermined boundary data on 8Dq to a volume integral of the 
inhomogeneous term j. Denoting the boundary data for u by u\q Do = f and n^(u) \sdo = fi, 
we have, explicitly, 

/ U/i - / rW)l = / (24) 

JaD L J "'-Do 

which are clearly necessary conditions on the boundary data for existence of a regular so- 
lution u. It is natural to ask whether such conditions are also sufficient. More precisely, 
assume that a continuous function j is given on D such that r 4 j is bounded at infinity. 
Give also two arbitrary continuous functions f and fi on dD which satisfy (24) for any 
choice of regular 7-harmonic function if>. We want to check whether there always exists a 
function u satisfying the Poisson equation A 7 w = j, with ru bounded at infinity and such 
that the boundary equations u\qd = fo, n^{u)\dD {) = fi are satisfied. The answer is yes as 
we prove next. Consider the Dirichlet problem A 7 w = j with u\qd = fo- Standard elliptic 
theory tells us that this problem always admits a unique solution u which tends to zero at 
infinity. Let us define fi on dD by f\ = n 7 («)|aD - Since u solves the Poisson equation, it 
follows from (23) that, for any regular 7-harmonic function ip, 

/ ipfi - fon*,(ip) dS 7 = / ipji]^. (25) 

JdD 1 J J D 

Our assumption is that f and fi satisfy (24) for any such ip. Subtracting (24) and (25) we 
JdD ^(f 1 ~ fi)^^ = 0- However, since we can take any regular 7-harmonic function if), 
the fact that the Dirichlet problem for the Laplace equation always admits a solution allows 
us to choose if)\dD to be any continuous function. This readily implies j\ = fi and hence 
that the overdetermined boundary data admits a decaying solution, as claimed. 

However, the compatibility condition (24) has a disadvantage, namely that it must be 
checked for an arbitrary decaying solution if) of the Laplace equation. This makes it useless 
in practical terms. Our aim is to reduce the number of solutions if) that must be checked in 
(24) in order to ensure compatibility of /o an d fi- Here is where axial symmetry plays an 
essential role. We restrict ourselves to axially symmetric functions fo and f\ on an axially 
symmetric boundary dD and we assume further that dD is simply connected, so that it is 
diffeomorphic to a 2-sphere. Using coordinates /i and (p on dD , let us denote by /i^d and 
fi jv (with us < A* at) the only two values of /i at which dD intersects the axis of symmetry. 
Let us also assume that they satisfy z (fis) < Zq{^n)) i-e. that /i increases from the "south" 
pole of the object (corresponding to fi,s) to the "north" pole (corresponding to /ijv). With 
this assumption, a direct calculation shows that dS^n^ = podfidtpn, where n is our usual 
normal vector n = —z d p + Pod z \oo- Similarly dS* 7 n 7 = e~ 4f/ °' E o p dfxdipn. For 7 = 7, (23) 
becomes, after a trivial angular integration, 

"Mat ^ r 

[if)n(u) -un(ip)}p Ldp = — / ipjr] y . (26) 

Z7T J Do 



I 
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For 7 = 7, (23) becomes 



rUN 
J us 



[ipn(u) — un(tl))]p e 4C7 °L dp, 



So 

US 

= ^zf ^'e- 12C V (27) 



2tt 



We want to choose a reduced set of functions ip for which the argument used above to show 
consistency of the data still holds. The following Lemma is probably known although we 
could not find an explicit reference for it. 



Lemma V.l Let h : [p Sl p N ] — > K. be a continuous function satisfying 

r " N h(p)dp 



\ yr ' r = 

Jus VPo(^) + ( Z M -y) 2 



for any constant y G (zg, z^), where z s = z (p s ) an d Zs = z (/1n). Then h = 0. 

Proof: Let us define a function ft, on <9-Do by extending h in an axially symmetric way, i.e. 
h(p, ip) = h(p). For any point q G E 3 let us define the function 

Uq) = [ jr^-^s,, 

JdD dist(-,g) 

where dist(-,g) denotes Euclidean distance between a point on dD and the point q. This 
function, also called "single layer potential", is well-defined throughout E 3 (including the 
axis of symmetry and the surface 8Dq [37]), is axially symmetric and vanishes at infinity. 
Moreover, potential theory tells us that Y h is C° everywhere (including dD ) and satisfies 
A 7 F fe = except on dD . The function h is directly related to the jumps of the first normal 
derivative of Yh on dD [37] (Yh can be physically interpreted as the potential created by a 
surface layer sitting on dD ). Thus h vanishes if and only if Yh is C 1 on E 3 . 

The hypothesis of the Lemma tells us that Yh vanishes on the piece of the axis lying 
between the south and the north pole of dD . Regularity of the Laplace equation shows 
that Y h is analytic except at dD . If Y h is identically zero in some neighbourhood of the 
axis inside dD , then continuity at dD and analyticity implies Yh = everywhere and 
h = would follow. Let us thus assume that there is a neighbourhood U of a point lying on 
the axis between the south and north poles in which Yh is not identically zero. Analyticity 
in Cartesian coordinates implies that Yh depends analytically on p 2 and z in cylindrical 
coordinates. Not being identically zero, there must exist a minimum value k G N and 
an analytic (not identically zero) function g(z) such that Y h \u = g{z)p 2k + 0(p 2k+2 ). The 
fact that Yh vanishes on the axis demands k > 1. Substitution into the Laplace equation 
dppY h + p~ x d p Y h + d zz Y h = yields Ak 2 g(z)p 2k ~ 2 + 0(p 2k ) = which is a contradiction. This 
completes the proof. □ 

This Lemma already suggests which subclass of regular 7-harmonic functions needs to 
considered for the compatibility of the boundary conditions. For the flat metric 7, a natural 
class of axially symmetric harmonic functions is 

z ) = , 9 , ) Tf> V e ( z s, z N ). (28) 
^p 2 + {z- y) 2 
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This expression is singular only at p = 0, z = y, which lies on the axis of symmetry outside 
Dq. Thus ip y is indeed a regular 7-harmonic function on Dq. 

For the 7 metric (corresponding to the Q-equations), we need to find a suitable class of 
axially symmetric solutions of A^u = 0. In order to find them more easily let us consider 
the semiplane IK = R + x H. defined as the subset {(f) = const., p > 0} of E 3 in cylindrical 
coordinates. Working directly on K has the advantage that axial symmetry is incorporated 
into the calculations from the very beginning. IK is endowed with a flat metric dp 2 + dz 2 . 
We choose the orientation so that *dp = —dz (and hence -kdz = dp). Obviously any axially 
symmetric function in E 3 immediately defines a function on K. Similarly a function on IK 
defines an axially symmetric function on E 3 . We shall use the same symbol to denote both 
functions (the precise meaning should be clear from the context). The field equations (15)- 
(16) can be translated into equations on K. It is straightforward to check that the 7-Laplace 
equation (i.e. the equation satisfied by Uq or Uq) becomes simply d(p*dUo) = outside the 
axis of symmetry. For the 7-Laplace equation, we first note the following simple identity, 
valid for any pair of functions f±, f 2 on IK, 

d/i A *df 2 = - * d/i A d/ 2 = - (d/i, d/ 2 ) 7 dp A dz. 

Thus the second equation in (15) can be rewritten, away from the axis of symmetry, as 

d (p * dQ' ) - 4pdn' A *dU = 0. (29) 

Our aim is to find suitable regular 7-harmonic functions, i.e. suitable solutions of this equa- 
tion on the domain IK E corresponding to the exterior domain Dq. More precisely, the surface 
S C E 3 projects into a line c in K defined parametrically as {z = z (p),p = p (p)}. This 
line separates IK into two regions, the exterior (denoted by K E ) and the interior. The as- 
sumption we made on the topology of S , namely that it has vanishing genus, implies that 
IK E is simply connected and that Co intersects p = at two values of z, namely z,s and zn- 
In order to determine suitable solutions of A^u = the following lemma is useful. 

Lemma V.2 For any y G (zs,z N ) and for any point (p,z) € K E define T y (p,z) G [0, 2tt) 
by 

cosT„(p, z) = — . , 

vV + (^-2/) 2 

sinT„(p, z) = — . % . 

Vp 2 + & - v) 2 

Then the PDE 

dZ y = cos T y dU + sin T y ★ dU , (30) 
with boundary condition lim p 2 +2 2^ 00 Z y = admits a unique solution on K E . 

Proof: The function T y is regular everywhere in K except at {z = y, p = 0} which lies 
outside K E . Similarly Uq is regular everywhere on this simply connected domain. So, in 
order to prove existence we only need to show that d(cos T y dU + sin T y -kdU ) = on K E . 
A simple calculation which uses only the fact that Uq is a flat-harmonic function on D Q gives 
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( i i A 

d(cos T y dUo + sin T y * dC/o) = — sin T y I dT y H — * dT y H — ) A dUo- 

Y tan 1 y p y 

Using the definition of T y above it is immediate to show that the term in parenthesis on the 
right hand side vanishes. Thus existence of Z y follows. Furthermore, asymptotic flatness 
implies that the right hand side of (30) tends to zero like 1/r 2 . Thus Z y is bounded at 
infinity and we can impose boundary data there. Since the general solution of (30) depends 
on an arbitrary additive constant, the lemma follows. □ 

Note that both T y and Uq are analytic functions on ~K E . It immediately follows that Z y 
is also analytic on this domain. Let us now define the function 



2U -2Z. : 



6 

z ) = 1 9 , ^ » y e ( z s, z N ). (3i) 

^ p 1 + {z- vy 

This function solves A^u = 0, as we show next. 

Lemma V.3 The function ^ y defined in (31) is a regular ^-harmonic function on D Q . 

Proof: We need to show (i) A^ y = on D Q , (ii) ty y admits a C 1 extension to dD and (iii) 
rty y is bounded at infinity. Property (iii) is immediate from the corresponding property of 
Z y . Furthermore Uq is C 1 on dD U D , which implies, from the PDE (30), that Z y is also 
C 1 on this subset. In addition is analytic except at the point {z = y, p = 0}; hence (ii) 
holds. In order to prove (i) it is sufficient to check the differential equation 

d(p* d^ y ) - Apd^y A *dU = 0, 

which is equivalent to = except on the symmetry axis p = ( A^^ = on the 

axis follows from the fact that ^ y is C 2 on D - in fact analytic). Note that ^ y e~ 2Uo+2Zy is 
just the flat-space harmonic function 1/r where r is the distance from p = 0, z = y. Using 
this and dZ y A -kdZ y = dUo A *dUo, a simple calculation gives the result. □ 

We can now state the necessary and sufficient conditions that the boundary data must 
satisfy so that A^u = j admits a decaying solution. 

Theorem V.l Let f , f 1 be continuous axially symmetric functions on a C 1 simply 
connected, axially symmetric surface S of E 3 . Let this surface be defined in cylindrical 
coordinates by {p = po(p), z = z (p), (f) — ip}, where p takes values in [p s , Pn] and ps < 
Pn are the only solutions of po{p) = 0. Call z$ = z(ps) and z^ = z(p^) and assume 
zs < z n (i-&- that these values correspond to the "south" and "north" poles of the surface, 
respectively). Denote by D the exterior region of this surface and let j be any axially 
symmetric function on D such that r A j is bounded at infinity. Let 7 be the flat metric 
and 7 = e~ 8(7o 7 7 where Uq is any regular flat-harmonic function on D . Then the Cauchy 
problem 

A7M = j, w| Eo = f , n(w)|s = h, 
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where n|s = —zod p + pod z , admits a regular solution if and only if the compatibility condi- 
tions 

7 = 7 : ft" bl> v fi ~ fo n(ip y )\ p | So dp = ± J Dq ^j^, G (z s , z N ) (32) 

7 = 7: [% fi - fo n(%)\ Poe-^ dp=l f Do Vyje-™^, Wy G (z s , z N )(33) 

are satisfied, where ip y and ^> y are given in (28) and (31) respectively. 

Proof: We give the proof for 7 = 7; the case 7 = 7 follows by setting Uq = everywhere. 
Necessity follows directly from Green's identity (23). In order to prove sufficiency, let u be 
the unique regular solution of the Dirichlet problem A^u = j, m|s = fo on D (which is 
known to exist). Define n(u)\z = f±. Green's identity and the fact that ty y is a regular 
7-harmonic function implies 

/ % h ~ fo n(* y ) Poe- 4U « d/i = — / ¥ w je- 12t V 
Jus 1 J s " 2ttJ Do 

Subtracting (33) we get J m ^ y (fi — fi)poe~ 4U ° dp = 0. We now apply Lemma V.l with 

~ Sq ^ 

the function h = e~ 2f/o-2Za | So p (A - /i)- It follows that h — and hence /1 = /1. Thus the 
Cauchy problem is solvable and the theorem follows. □ 

For the first order perturbations, the inhomogeneous term j vanishes and this theorem 
provides necessary and sufficient conditions involving the boundary data only (and hence 
conditions on the interior perturbations via the perturbed matching conditions described in 
the previous section). 

For the second order perturbations things are not so easy because the equations are 
inhomogeneous (j 7^ 0). In principle, one would need to integrate the first order functions 
Uq, Q' in order to compute j and thence J Do ipyji]^ (and the corresponding expression for 7). 
In some practical situations, an alternative procedure would be finding a particular solution 



Up of 



A^u = j. 



so that the homogeneous compatibility conditions (i.e. with j = 0) can be applied to the 
function Uh = u — u p , which solves the homogeneous equation A^Uh = 0. More specifically 
we should need to check whether the Cauchy data {/ — u p \^ Q , f\ — n(u p )\-£ } (where / 
and fi are the Cauchy data for u) satisfy the homogeneous compatibility conditions in 
theorem V.l. Obviously this approach relies heavily on knowing a particular solution of the 
inhomogeneous equation. 

It is clear that, generically, we shall not be able to integrate the first order equations 
explicitly or find a particular solution of the inhomogeneous equation. We should still like 
to be able to treat the problem in a satisfactory way. The key idea which enables us to do so 
is to rewrite the volume integrals on the right-hand sides of (32) or (33) as surface integrals, 
and hence rewrite everything in terms of boundary data (or integrals thereof). Suppose that 
we were able to find an axially symmetric vector T = T p d p + T z d z on Do such that (with V a 
denoting covariant derivative with respect to 7) V a T a = ip y j (for 7) or V a T a = ^ y je~ 12Uo 
(for 7). Then we could use Gauss' identity to transform the volume integral into a surface 
integral on E . We now show that this is indeed possible. As before, it is useful to work on 
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the two-dimensional space K. E . Denning the one-form T = T p dp + T z dz we can translate 
the equation above for T into an equation for T on K. E by means of the general identity 
d(p * T) = -pV a T a dp A dz. 

We start with the equation for Uq in (16). We obviously have j = — e~ 4U °(dQ' , df2()) 7 so 
that the equation for T reads, using the explicit form (28) for ip y , 

d (p * T) + sin T y e^ Uo dn' A *dfi(, = 0. (34) 

Our aim is to find a solution of this equation. In order to do this, the following Lemma 
turns out to be useful. 

Lemma V.4 Let Q' Q be a solution of the second equation in (15) and Z y be defined as in 
(30). Then the equations 



dS 1 = e - 2Uo+2Z y [- (1 + cos T„) d% - sin * dQ' ] , 
dS 2 = e - 2U °- 2Z y [(1 - cos T„) dtt' - sin T y * dtt' } , 
admit unique solutions which tend to zero at infinity. 



(35) 



Proof: The solutions, if they exist, are unique except for additive constants. Asymptotic 
flatness of Q' implies that Si and S 2 each tend to a constant at infinity; these additive 
constants can be chosen so that Si and S 2 vanish at infinity. Thus the only non-trivial 
part of the proof is to show existence of the solutions. Simple-connectedness of K B implies 
that we only need to check oj = d ( e - 2(7 o+25z y _^ cog f y ^Q' Q _ s j n f y * dtt' ]) = 0, where 
5 = ±1. A straightforward, if somewhat long, calculation using the equations for Uq, Q' 
and Z y gives 



w = e MZ »- 2t/0 siiiT 1/ 



^dT,, dz 
dT y + —^r + — 
tan l y p 



A d% = 0, 



where again the explicit expression for T y is used in the last equality. □ 

The equations for Si and S 2 imply dSi A dS 2 = 2 sin T y e~ 4U °dQ' A -kd£l' , which is the 
crucial fact allowing us to solve (34). Indeed, defining 

Tl = ^Si*dS 2 , (36) 
2p 

it is immediate to check that (34) is satisfied for T = T x . This expression is apparently 
singular at p = 0. However, on the axis of symmetry we have cosT^ = +1 or cosT^ = — 1 
depending on whether we are above the north pole or below the south pole, respectively. 
Consequently, we have Si = on the subset of the symmetry axis below the south pole 
(from the equation it satisfies and the fact that Si vanishes at infinity) and dS2 vanishes 
on the axis above the north pole. Combining this with analyticity of Q' everywhere on D 
(including the axis), regularity of T\ follows. 

Considering the equation for f2' ', in this case we have j = 8 (dQ' , dU^)-. Using the explicit 
expression (31) for ty y , the equation we need to satisfy is 

V a T a = - 8 e -iUo-2Z y W) du , } 

^p 2 + {z-y) 2 
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which, in terms of exterior forms, reads 

d(p-kT) =8 sin r y e- 2U °- 2Z «dn' A *dU' Q . (37) 

We now use the fact that U$ is a flat-harmonic function. Thus (compare Lemma V.2) we 
can uniquely define a function Z y on K E by the equation 

dZ' = cos T y dU' Q + sin T y * dU' , lim Z' = 0. 

It is then straightforward to check that 

T 2 = --S 2 *d(Z' y + Ufi (38) 
p 

satisfies (37). Again, regularity at the axis follows because S2 vanishes on the axis above 
the north pole and d(Z y + U' ) is zero on the axis below the south pole. 

We end this section by summarizing its main results in the form of the following Theorem. 

Theorem V.2 Let the assumptions and notation of theorem V.l hold. Then 

(i) the Cauchy boundary value problem 

A 7 ^ = 0, ^k = /o, n(£/o)k = /i, 
admits a regular solution on D if and only if 

/ bPy fi - fo n(4>y)} Po| Eo dp = 0, Wy E (z s , z N ), 

(ii) the Cauchy boundary value problem 

A,n' - 4 (dtt' , dV Q \ = 0, Q' \x = fo, n (%) | So = f u 
admits a regular solution on D if and only if 

[*„/i " /o^(^)]poe- 4% | So dp = 0, \/y E (z s , z N ), 

(Hi) the Cauchy boundary value problem 

A 7 ^' + e- 4l/ ° (dn' , dtt'^ = 0, Etf| Eo = fo, n(U%) |s = /1, 
admits a regular solution on D if and only if 

'•Mat 

[ipy fi ~ fo n{ip y ) - n (Ti)] po| So dp = 0, \/y E (z s , z N ), 

'MS 

and (iv) the Cauchy boundary value problem 

AX - 8 (&V , dU' Q \ - 4 W, dU \ = 0, ^'| So = f , n (fig) | So = f u 
admits a regular solution on D if and only if 

/ h - fo n(* y )) e~ 4Uo - n (T 2 )] p | So dp = 0, My E (z s , z N ), 

•J 

where ip y , ^ y , 2~i and T 2 are given in (28), (31), (36) and (38) respectively. 



f 

J Uf- 
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Remark. Since writing down Ti and T2 requires the integration of d^ it may seem at 
first sight that one has not really gained anything with respect to the volume integral in 
theorem V.f . The difference however is substantial because we only need to know S2 on the 
boundary S . So by projecting equation (35) into the boundary we get an ODE for S^ls,, 
which can in principle be solved by quadratures. Hence the problem reduces to perform- 
ing integrals, which is of course much easier than solving PDEs. Thus the compatibility 
conditions are truly written in terms of the boundary data alone, as we wished. 



VI. PERTURBATIONS AROUND SPHERICALLY SYMMETRIC STATIC 
BACKGROUND CONFIGURATIONS 

Everything we have discussed so far holds for any stationary and axially symmetric per- 
turbation of a static and axially symmetric background. Of course in many cases of physical 
interest, and especially for perfect fluids, one expects that equilibrium (non-rotating) con- 
figurations of isolated bodies are spherically symmetric. In the case of perfect fluids this 
has been rigorously proven for a large class of equations of state [12]. It is therefore of 
interest to specialize the previous results to the case when the background spacetime is in 
fact spherically symmetric. This implies in particular that the exterior background metric, 
being vacuum, corresponds to the Schwarzschild metric 



ds 2 = - (l - ^ dt 2 + - ^ 1 dr 2 + r 2 {d9 2 + sin 2 9 d0 2 ) , 

where m is the total mass of the spacetime. The static Killing vector is given by £ B = d t , 
and we fix the axial Killing vector to be fj^ = in these coordinates. (The non-uniqueness 
in the choice of axial symmetry in a spherically symmetric situation is precisely what allows 
us to orient the coordinate system so that this vector is the limit of the axial symmetries of 
the perturbed spaces.) From the definition of U we get 

t/i^i-^y (39) 



2 V r J 

The intrinsic definition of the Weyl coordinate p, p 2 = -(£f ',£<?) g E(f% \fjg) g E + (£f ,Vo) 2 g E, 



gives 



(40) 

In order to determine z we use (dz,dp) g E = 0, (dz, dz) g E = (dp, dp) g E, which can be solved 
to give z = z ± (r — m)cos9, where zq is a constant. With these coordinate changes, 
{r, 9, 0} can be regarded as coordinates on (E 3 , 7). Choosing {dr, d#} to be positively 
oriented on a plane = const., i.e. that *dr oc d# with a positive proportionality factor, 
and choosing (as we have above) that {dz, dp} is also positively oriented, the + sign above 
is selected. We can further choose z = without loss of generality, so we have 

z = (r — m) cos 9. (41) 
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Spherical symmetry of the whole background spacetime requires that the surface of the 
body Xjf is defined by r = ro(> 2m). The embedding for this surface can be chosen to be 
Xq : {t, ip, /i} — > {t = r, (ft = ip, r = r , 9 = n — //}, with G [0, 27r) and /i G [0, 7r]. This 
choice for /i is motivated by the fact that it increases from the south to the north pole of 
the body, as required. In fact, we have /is — and /xjv = vr. The range for the constant 
7/ introduced in the previous section is then given by — r + m < y < r Q — m. With this 
embedding we clearly have 



/ 2m 

p (fi) = r sin/iWl , z (/i) = -(r - m) cos// (42) 

V ?"o 

The vector e = dx^(<9 M ) has norm r 2 , which implies that n (which is defined to have the 
same norm) reads 



no = -r \/l — d r 

r 



(43) 

So 



We can now write down the first and second order perturbed Ernst equations in Schwarz- 
schild coordinates. Performing the coordinate transformation (40), (41) the flat metric 7 in 
E 3 becomes 



7 = (r ' m) ^ C ° S 9 + rW ) + r (r - 2m) sin 2 M0 2 , 



r 



which implies the following form for the first order Ernst equations (15) 

cos 

r (r - 2m) % )TT + 2(r - m)U' Q , r + C/^ fl + — —U' Qfi = 0, 

sin t/ 

cos 

r (r - 2m) fi' , rr + 2(r - 3m)fi' , r + + — — % fi = 0. (44) 

sin u 

The second order equations (16) become 

n Jn'2 „2n/2 

r (r - 2m) U , rr + 2(r - m)U , r + U fie + —U , e + + - _ = 0, 

cos 

r (r - 2m) (Q' \ rr - 8fi' , r ^, r ) + 2(r - 3m) ^', r + fi' > + — - 85^1^ = 0. 

Sill u 

Let us analyze the equation for U' Q . Expanding in axisymmetric spherical harmonics 

00 

0S = $>, (1) (r)fl(cos0), 

we get a collection of ordinary differential equations which, after performing the change of 
variables r = m(l + x), read 

(l-x 2 )«£ B -2 a « I ( 2 + Z(Z + l)« I (1) = 0. 
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These equations can be solved as u~p = aiP t (x) + biP t (x) f*^ nz^fp^j ? where Pi(x) is the 
Z-th Legendre polynomial and a\ and q are constants (note that the integral in the second 
summand converges). Asymptotic flatness demands that U' Q tends to zero when x — > oo. 
Thus ai = for all / and the most general solution for U' Q is 

P^ = f>P(cos0)p(--l) / m 1 ^L^, 

^ ;y +oc (l-y 2 )P/(y) 

This expression also gives the general solution of the homogeneous part of Uq. For Q' Q , the 
expansion in spherical harmonics Q' = Ym^o w i 1 \ t )^( cos ^) transforms (44) into 

(1 " * 2 ) w£L - 2 (x - 2) «,£ + /(/ + !) = 



where again the x variable has been used. The general solution is now 

wjp = qP z (_2 ' 2) (x) + d,.P/~ 2,2) (x) J? - — fe^Sr: , where P/ _2 ' 2) (x) is the Jacobi polyno- 

mial and q and d\ are constants (we are using the notation of [38]). Imposing the condition 
that the solution tends to zero at infinity we conclude that q = and we can write the 
general solution as 

r \ 1 (y — ^)dy 



^ " Vm d/ + i) 3 W ( " 2,2) (y)) 2 

Having obtained these solutions, we consider next the boundary data that these functions 
must satisfy. In Proposition IV. 1 we obtained expressions that involve only tangential deriva- 
tives of several scalar objects. We need to evaluate them. First of all we notice that spherical 
symmetry implies that 

e(P |s ) = e(n(Po)|s )=0, (45) 

which substantially simplifies the boundary conditions. The explicit forms of n(Uo)\s , X 
and X\ are directly obtained from their definitions as 



n(U )\ 



-1 



So 



x ^/x[ 



„ v .. „ - cos/ism/i 

— Q 5 ^1 



Xq — cos 2 [i x\ — cos 2 \i ' 



where x is the value of x on the surface, i.e. x = r /m — 1. Since the background spacetime 
satisfies the matching conditions, we obviously must have Vq = l/21og((rro — l)/( x o + 1)) 



and nVo = — \j\/x\ — 1. The functions Pi, Qi, P2 and Q2 are still arbitrary because they 
determine how the unperturbed surface is deformed to first and second order. The same 
is true regarding the functions Vq, HVq, V ", nV " for the P-boundary conditions and Wq, 
uWq, Wq, HWq for the 0-boundary conditions, which depend on the interior perturbations. 
Using all these expressions, the boundary conditions described in Proposition IV. 1 simplify, 
for spherical backgrounds, to 
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OlSo 



n(^)|E =n\/ ' + 



dP 1 

dfj, 



a; 



Vq" -2P 1 -^-2Q 1 nV^- ( 2 



d(PiQi) 
d\x 



-Q2 + QI 



y/^l-l) \/x\ - 1 ' 

rr 



(Pf - Q?) ^oa/^o - 1 + 2Pi<5i cos /x sin /x 



xl — COS 2 /X 



SWlla, = W„" + 2- 



d/x \ <i/x 



_ 2 d(j^avg) + 



(i/x 



cos /x c?V^' 
sin 



x 



in/x d/x - 1 



p , (Pi 2 -Qi)cos/xsin/x rr 



— COS 2 /X 



1 d / 2PiQia;o v a^o — 1 si 11 / 1 



sin /x dfx 



Xq — cos 2 /x 



+ Q\ cos /x) 



^0 



Xq — COS 2 /X 



sin /x y — 1 



W^, n(fi' )| So =n^, ^| So = ^'_ 2 P 1 ^-0-2g 1 n^, 



a/x 



n ^0 
wi —, — 

o/x 



a/x 



cos /x cWq a;o — 4 



sin/x d/x a/^^ - 1 



Our next aim is to find what compatibility conditions these data must satisfy in order to 
admit asymptotically flat solutions. In particular we need to find the functions ip y and ^f y 
defined in the previous section. Recalling the definition ip y — 1/ p 2 + (z — y) 2 and using 
(40) and (41) we get the explicit form 



a/ m 2 x 2 + y 2 — 2mxy cos 9 — m 2 sin 2 9 



where y is a constant satisfying \y\ < mx . For *ff y we first need to integrate the partial 
differential equation (30) for Z y . Contracting this equation with d$ and using the fact that 
Uq depends only on r we get 



msin9 



96 a/ m 2 x 2 + y 2 — 2mxy cos 9 — m? sin 2 9 
which, together with the boundary condition at infinity, gives 



e -z y 



Zp(x) 
(y — m) \J x 2 — 1 



yx — m cos 9 — y m 2 x 2 + y 2 — 2mxy cos 9 — m 2 sin 2 9 
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where Zq is an arbitrary function of x obeying Zq — > 1 as x — > oo. (Using L'Hopital's rule 
one can easily check that this form for e~ Zy has the correct limit at the special value y — m.) 
Then contracting (30) with d r we find that Z is constant. Thus 



e -z y = 



yx — m cos 9 — yj m 2 x 2 + y 2 — 2mxy cos 9 — w? sin 2 9 



(y — m) \J x 2 — 1 

We can finally write down the explicit expression for ty y , which reads 



^yx — m cos — yj m 2 x 2 + y 2 — 2mxy cos 9 — m 2 sin 2 
(y — m) 2 (x + l) 2 ^ rn?x 2 + y 2 — 2mxy cos 6> — m? sin 2 



Having obtained ^ and ^ explicitly, the compatibility conditions that the boundary data 
for Uq, Q' , Uq and Qq must satisfy in the spherically symmetric case are direct consequences 
of theorem V.2 in the previous section. 



Acknowledgments 

This research was primarily funded by EPSRC grant GR/R 53685, which supported R. 
Vera as a Research Assistant and M. Mars as a Visiting Fellow. We are grateful to the other 
Visiting Fellow supported by the grant, Prof. J.M.M. Senovilla, for helpful discussions. Dur- 
ing the final stages of preparation, R. Vera was supported by the Irish Research Council for 
Science, Engineering and Technology postdoctoral fellowship PD/2002/108, and later by a 
Basque Government grant Ref. BFI05.335. M. Mars wishes to thank Robert Bartnik for 
pointing out the connection of the problem with the Dirichlet-Neumann map and for provid- 
ing useful references. He also acknowledges financial support from the Spanish Government 
under Project BFME2003-02121. Many of the calculations were verified with the help of 
the computer algebra system REDUCE. 



APPENDIX: THE PERTURBED BOUNDARY CONDITIONS 

In section IV we introduced four functions V € , nV t , W e and nW e such that the full matching 
conditions for the Ernst potential are U e \-£ e = V e , n t (U e )\Y, t = nV e , fi e |s e = W e , n e (f2 e )|s e = 
nW e . We want to take first and second derivatives with respect to e (at fixed //) in order to 
obtain the perturbed matching conditions in Proposition IV. 1. Since the calculations for U € 
and Q t are very similar, let us introduce a function F e (p, z, e) which satisfies 

F e \z t = H e , n e (F e )| Se = J e , (A.l) 

for some functions H e (fi, e), J £ (/i, e). We shall obtain the perturbed expressions for this 
generic function and then specialize to U e and to Q t . To do so we need to differentiate 
(A.l) with respect to e. The right hand sides are defined as functions of (//, e) so that for 
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them d/de = d/de, but on the left hand sides we have functions such as F e (p, z, e)|s e = 
F e (p(fi, e), z(fi, e), e). For the latter we have to take a convective derivative with velocity 

— * — * 

Z 1>e , defined by the obvious generalization of (18) when we evaluate at S e instead of S . Z 2tt 
is defined similarly. 

Evaluating (A.l) at e = 0, we immediately get the unperturbed boundary data Po|s = 
H , n(F )| So = Jo, and taking the first convective derivative of P e gives (on simple rear- 
rangement) 

F^ = H' - Pi^FqIeo) - Qin(F )U . 
To obtain the first derivative of n e (F e ) = n\diF € we calculate 



Ik 
dn\ 
~de~ 



dF 



diF e + Z\ e nididjF e + n%F' t 



Substituting for Z\ and from (18) and (19), using the identity 



(A.2) 



which like subsequent similar identities follows from integration by parts using <9 M (e) 
Xie — X n and d^{n) = X e + A^n, and rearranging, we obtain, at e = 0, 



n(Fo)k = Jo-A^(W|Eo)-Qi(nV9 i a,-Fo| E o) 



(A.3) 



+ 



The second derivative of F e follows by applying the convective derivative twice. Thus we 
obtain 

K = {Fl + ziAF,) 1 + ziMK + ziAjFe) 

= F'l + 2ZIAF* + ZlAFe + Z\ I Z{J);0 J F I . 



— * — * 

Substituting the values of Z\ and Z 2 from (18), using the identities (A.2) and 

eVc^F | So = ^(F |eo) - (^i^(F | So ) - X n(F )| So ), 
and rearranging we obtain 



(A.4) 



Fn 



I So 



H'> - 2P 1 d fl (F^ - 2Q 1 n(Fo)ko ~ Qi (nWdftF^) 
-2P 1 g 1 9 M (n(Fo)| Eo ) - ^^(FqIeJ 

- (Q 2 + P*X Q - 2P 1 Q 1 X 1 ) n(P )k - (P 2 - P^ - 2PiQiX ) ^(P | E() ), 
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The last of these expressions that we need is the second derivative of n e (F t ). This we can 
obtain by calculating 

[*W)He. = ^0^ + 2^(8^' + Zi e d i d j F e ) 

+nl(ZlAd 3 F e + Z{ e Zldi d AFe + ^[^AK + W)| E< 
The terms in F" and F' t on the right at e = are readily evaluated on substituting for ^ 

— * 

and Zi from (19) and (18), and again using the identities (A. 2) and (A. 4). They lead to 
r1( /•;;'):>:„ and 

2Q 1 (nWdidjF^) + 29 M (P 1 n(^)| Eo ) - (2^- + 2Q 1 X 1 ^j d,(F^ + 2Q 1 X n(F^ 

respectively. To evaluate n l Z{ZididjdkF \-£ in the form we require, we will need the iden- 
tities 

e i e i n k d i d j d k F \z = «9^(n(F )| So ) - X <9 M (F | So ) - X 1 d ll {fi{F Q ) | Eo ) - X ^(F | So ) (A.5) 

-Xin(F„)| So - XoeVM-Fota - 2X 1 eWd i a,-F | Eo + X nV,9,^F | So , 

eWri^c^Fok = <9>W<9^F | So ) - 2X d^n(F )\^ ) - 2X 1 (nV9 i a,-F | Eo ) (A.6) 

+2X (X ^(F | So ) +X 1 n(F )| Eo ). 

Using these, the previous identities (A. 2) and (A. 4), and the values of Z±, Z 2 , and 
from (18) and (19), we obtain 

^TOIeo = ^-2Q 1 (n i nW i d j F^ ) - 2^(P 1 n(F^)| Eo ) 

2^ + 2Q!X^ ^(^)| So - 2Q 1 X n(F ( ;)| Eo 
-Q? (nVn^a^ok) - 2P 1 Q 1 d, (nV<9^F | So ) ~ P?d^(n( F o)ko) 
- (Q2 + 2g?X + 2Q^ + P^Xq - 2P 1 g 1 X 1 ^) (nV<9^Fo| S() ) (A.7) 

+ (-P 2 + P2l ) + + 2QlX^j d,(n(Fo)U 

+ ( 2p i^ + 2PiQiX^ «9^(P | So ) 



+ 



dP 2 + gpj%) _ ^d(Qi) + 2 Pl x rfgi 



+ 



d/j, d\x d\x d\i 

-Q 2 X - P?X 2 + 2P 1 Q 1 X X 1 - 2Q\X\] n(Po)|s 

% + d(p*x ) _ x d{Q\) _ 2PiXi dQi 

d\x dji djj, djj, 

+Q 2 X 1 + P*X Q X X - 2P\Q\X\ - 2Q 2 1 X X 1 ] <9 M (F | So ). 

These expressions are in fact identities which hold for any axially symmetric function, in 
particular for U e and Q e , but they are not in a form given by the boundary data, since they 
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involve {n l v? 'didjFo\-£ ) and {n l n^n k didjdkFQ\Y,^) ■ However, because the functions Uq, Uq, fl' Q 
satisfy the elliptic equations (14) and (15) which relate second order tangential to second 
order normal derivatives we can eliminate the second and third order transverse derivatives. 
Since 7%o = ^("V + eV) + ^(fy)* W|s , So{d pp + d zz )F = (nW + eV)d^F and we 
immediately find the identity 



(nV^Folso) =5o(A 7 Fo)| So -^(Fo|E )+X 1 ^(F | s J-X n(F )| E( -^ ( 9 M (Fo| s J + ^n(F )| So . 

Po Po 

(A.8) 

Applying this identity to Uq, U' and Q' and using the field equations they satisfy we obtain 
the following expressions for the transverse-transverse derivatives on the boundary 

(wW^E/oIe,,) = -^(E/qIeo) + - ^ ^(C/o|so) - (a - n(C/ )|s , (A.9) 

(nW$0,-^| Eo ) = -^(USIeo) + (*i - ^ 9,(^| So ) - (x - ^ n(^)| So (A.10) 
(nV^Q | So ) = 4nTO| So ^o)|Eo+4W| So )^(C/ |Eo)-^(Q( ) | So ) + 

+ - ^ ^(fi | So ) - (x - ^ n(^)| So . (A.ll) 

It only remains to evaluate (n l n 3 n k didjdkUo\j] ) in terms of known boundary data. In order 
to do that we need a similar identity but now involving third derivatives. A straightforward 
calculation gives 

(n i n*n k d i d j d k F \x ) = S n (A 7 F ) | So - «9^(n(F )| So ) + 2X 3 wt (F | Eo ) 

+ (3X, - ^ ^(n(F )| Eo ) + (-X + j^j (nV^^Folso) 

+ ( A — 3X Xi + — X — ) ^(-^o|eo) 

V Po Po / 

+ (x x - 2X1 + Xl + ^ + 4) n(F )| Eo , 

V Po p J 

which, again, holds for any axially symmetric function in E 3 . Applying it to Uq, and using 
AjUq = and (A.9), we obtain 

(rjnWdidjdkUol^) = -^(n(C/o)|Eo)+9 w (C/o|E ) (sX - + 

+WIe ) fx - 4* *i + 2^X + ^ - 2 -^p) 
V Po Po Po / 



+n(C/ )|E (*i - 2X 2 + 2X 2 + ^X, - ^X + (A.12) 
V Po Po Po / 
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Substituting the transverse-transverse derivatives (A. 9), (A. 10), and (A. 12) into (A. 3), (A. 5) 
and (A. 7) with the substitutions F — > U, H — > V and J — > nV we find the Cauchy 
[/-boundary data in Proposition IV. 1. Note that in this evaluation we need to use the 
formulae for U' q \y, and n(U'^)\-z in those for Uq\^ and n(C/o)|s - Finally, the substitution 
F — > Q, if — > and J — > nIV and use of (A. 11) and the fact that f2 = gives us the 
0-boundary data. This completes the proof of the Proposition. 
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